Polarization dependence of spin excitations in BaCu 2 Si207. 
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The polarization dependence of magnetic excitations in the quasi one-dimensional antiferromag- 
net BaCu2Si207 is studied as a function of momentum and energy transfer. The results of inelastic 
neutron scattering measurements are directly compared to semi-analytical calculations based on the 
chain-Mean Field and Random Phase approximations. A quantitative agreement between theoreti- 
cally calculated and experimentally measured dynamic structure factors of transverse spin fluctua- 
tions is obtained. In contrast, substantial discrepancies are found for longitudinal polarization. This 
behavior is attributed to intrinsic limitations of the RPA that ignores correlation effects. 



I. INTRODUCTION 



Excitations in weakly ordered quasi-one-dimensional 
(quasi-lD) antiferromagnets (AFs) are a topic of consid- 
erable current interest in the field of quantum magnetism. 
Particularly intriguing is the problem of the so-called lon- 
gitudinal mode (LM) , a magnon excitation polarized par- 
allel to the direction of ordered moment. The discovery 
of a coherent LM in KCuF3p(Refs. |l|,^|) confirmed pre- 
vious-theoretical predictions^ based on the chain-Mean 
Fieldo (chain-ME-i and Random Phase Approximation 
(RPA) theories m Currently chain-MF/RPA indeed ap- 
pears to be the most versatile analytical framework for 
treating weakly-coupled quantum spin chains. However, 
the KCUF3 experiments also highlighted certain limita- 
tions of this approach. In particular, the chain-MF/RPA 
can not, by its very definition, account for the experimen- 
tally observed finite lifetime (broadening) of the LM. 

In a recent short paperi we reported polarization- 
sensitive neutron scattering measurements of the dy- 
namic spin structure factor in another model quasi- ID 
antifcrromagnet, namely BaCu2Si20y. This S = 1/2 
system has much weaker inter-chain interactions and 
low-temperature ordered moment than KCUF3. Pre- 
liminary results indicated that, unlike in KCUF3, in 
BaCu2Si207 there is no well-defined longitudinal mode. 
Instead, the longitudinal spectrum is best described as a 
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single broad asymmetric continuum feature. This stark 
discrepancy with the predictions of the chain-MF/RPA 
model came as surprise. Indeed, for the transverse- 
polarized spectrum of BaCu 2 Si20y, earlier neutron scat- 
tering work confirmed excellent agreement with chain- 
MF/RPA theor« n at least as far as excitation energies 
were concerned.lZrO The apparent paradox is not fully re- 
solved to date. This is in part due to that only very 
limited data are available for longitudinal-polarized ex- 
citations. Even for the transverse-polarized spectrum, 
the existing wealth of high-resolution neutron data could 
not be quantitatively compared to theoretical predictions, 
for lack of calculations based on the specific geometry 
of inter-chain interactions in BaCu2Si207. The present 
work addresses both these issues and involves a detailed 
experimental and theoretical study of the polarization 
dependence of magnetic excitations in this compound. 
First, we further exploit the technique of polarization 
analysis described in Ref. H to investigate the wave vec- 
tor dependence of longitudinal excitations. We then per- 
form chain-MF/RPA calculations of the dynamic struc- 
ture factor for the exchange topology and constants of 
BaCu2Si2C>7. This enables us to perform a direct quan- 
titative comparison between theory and experiment for 
both energies and intensities of the coherent and diffuse 
components of the dynamic spin correlation functions. 

Magnetic interactions in BaCu2Si207 have_heen previ- 
ously thoroughly studied using bulk methods Jj'E3.aeutron 
diffractionJj'El and inelastic neutron scattering£loEj'E3 

The silicate BaCu2Si207crystallizes in an orthorhom- 
bic structure (space group Pnma, a = 6.862 A, b = 
13.178 A, c = 6.897 A) with slightly zigzag AF 5 = 1/2 
chains of Cu 2+ ions running along the c axis. The in- 



FIG. 1: Typical constant-_E scans (a,b) and constant-Q scans 
(c) measured in BaCu2Si2C>7 in magnetic fields H — 1.5 T 
(open circles) and H = 2.2 T (solid circles) applied along the 
crystallographic c-axis. The dashed lines in (a) and (b) repre- 
sent the background obtained by linear interpolation between 
intensities measured at I = 0.8 and I = 1.2. In (c) the dashed 
line is the background scan measured at q = (0, —0.5, 1.2). 

respectively.00 The transition involves a re-orientatiaa 
of the ordered staggered magnetization in the system.Eil 
As explained in Ref. ^, this leads to a drastic change in 
the polarization-dependent part of the scattering cross 
section for unpolarized neutrons. The effect on the scat- 
tering intensity from longitudinal (parallel to the ordered 
moment) and transverse (perpendicular to the ordered 
moment) spin fluctuations is different, which allows us 
to separate the two components. In general, the mea- 
sured intensity can be expanded as: 

I(q,w) oc S'- L (q,w)(l + cos 2 a q ) 

+ 5 ll (q,w)sin 2 a q + B(q,u/). (1) 

In this equation S ± (q, w) and S"(q, to) are the magnetic 
dynamic structure factors for transverse and longitudinal 
polarizations, respectively. The wave vector dependent 
angle a q is measured between the momentum transfer 
q and the direction of ordered moment. The orienta- 
tions of the latter was previously determined using neu- 



FIG. 3: Transverse (a) and longitudinal (b) components of a 
constant-Q scan measured in BaCu2Si2C>7 at q = (0, —0.5, 1). 
Lines and the plot shown in (c) are as in Fig. bl 

detail in Ref. ||. 



III. EXPERIMENTAL RESULTS 

Typical raw data sets measured in constant-Q and 
constant-i? modes at T = 1.5 K are shown in Fig. [|. 
At energies in excess of about 2 A the scattering is prac- 
tically unaffected by the phase transition (Fig. |l]b). The 
contrast in inelastic intensity measured at two different 
field values is most apparent at energy transfers of about 
A( Fig. [l|a and c). Separating the longitudinal and trans- 
verse contributions as described in the previous section 
yields the constant-Q scans shown in Figs. || [I| The evo- 
lution of the instrumental FWHM resolution ellipsoid in 
the course of each scan is shown in the right part of each 
figure. Typical constant-E 1 data are shown in Fig. ||. A 
contour and false color plot based on a series of 10 such 
scans taken with 1 meV energy step is shown in Fig. ^. 

Certain important features of the measured transverse 
and longitudinal dynamic structure factors can be iden- 
tified even without a quantitative data analysis. An im- 
portant experimental observation is that longitudinal ex- 
citations show a steep dispersion along the chains. As 
can be seen in Fig. |j, the corresponding spin velocity 
is the same as for transverse-polarized spin waves. Fur- 
thermore, at high energy transfers (above « 7 meV) the 
scattering is almost polarization-independent to within 
experimental accuracy and resolution (Figs. ||-|5|). Such 
behavior is consistent with our general expectation that 
inter-chain interactions become almost irrelevant at ener- 
gies well above the gap energy A. The dynamic structure 
factor in this regime is as in isolated chains, and is there- 
fore almost isotropic. 

At smaller energy transfers the structure factors for 
longitudinal and transverse polarizations are noticeably 
different. As observed in previous detailed studies,El 
transverse-polarization constant-Q scans are character- 



FIG. 6: Contour and false-color plot of the transverse- 
polarized (left) and longitudinal-polarized (right) inelastic 
scattering measured in BaCu2Si207 near the ID AF zone- 
center (0, 0, 1). 



ized by a sharp spin wave peak, whose position and 
intensity is strongly dependent on momentum transfer 
q± in the direction perpendicular to the spin chains. 
The effect of this pronounced dispersion can be seen 
in Figs. 2a-4a. fn contrast, longitudinal-polarized scans 
lack the sharp component and are almost independent of 
q± (Figs. 2b-4b). Such behavior is reminiscent of that 
for the transverse-polarized continuum that also shows 
very little variation with qj_ . 8 



IV. THEORY 



Before discussing the quantitative analysis of the ex- 
perimental data we shall describe the application of the 
Urn (meV) chain-MF/RPA approach to the problem of weakly cou- 




I (r. I. Ll) The Fourier transform of the inter-chain coupling is de- 

fined as 

J'(q) = J x cos(q x ) + J y cos (q y ) 

+ J 3 [cos(q x + q y ) + cos(q x - q y )\ . (3) 



i,j n 

h = 2(J X - J y - 2J 3 ) rrio = J f m . (7) 

The Hamiltonian (Q) describes an ensemble of uncoupled 
spin-^ Heisenberg chains in a staggered magnetic field 

ffid = ^JS n -s„ + M-i) n s z n+l . (8) 



n 



The next step is to find a solution for an isolated chain in 
an external field h. Since in the limit of weak inter-chain 
coupling the latter is expected to be small compared to 
J, it is possible to determine dynamical correlation func- 
tions at low energies Tiuo <C J by means of field theory 
methods. A standard bosonization analysis gives the fol- 
lowing scaling limit of (I 



7~^id = I dx 



v 



~ {dx<t>) + 7T + Chcos{V2i<t>) 



(9) 



In this formula v = TiJap-12 is the spin velocity of the 
spin- 1/2 Heisenberg chainE3 and C is a non-universal con- 
stant that was calculated in Ref. [lj. The model (||) is 
known as the quantum Sine Gordon model (SGM) and 
is exactly solvable. The spectrum is formed by scatter- 
ing states of four particles, called soliton s, antisoliton s, 
breather B\ and breather B^. The breathers themselves 
are soliton- antisoliton bound states. All foou: particles 
have gapped relativist ic dispersion relations :E3 

A 

E a = A cosh 6 , P (y = — sinh 9 , a = s, s, B\ , 

v 

E B2 = V^A cosh (9, Pb 2 = ~ sinh (9. (10) 

v 
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FIG. 8: Comparison of simulated scans across transverse- 
polarized continuum at q = (0, 0, 1) based on the exact chain- 
MF/RPA result (solid line) and the empirical TMA fitting 
function that was used to analyze the neutron scattering data 
(dotted line). Given the effects of experimental resolution 
that were taken into account in these simulations, the two 
curves are virtually identical. 



Here (q) is the spin wave dispersion relation given by: 

wi(q) = y j2 sin V) + jJTj [I J 'l + 2J/ (1)] > ( 24 ) 

where J'(q) is defined by Eq. |[ 

The second component of the fit function for trans- 
verse excitations approximates the continuum. We have 
previously found that, at least for wave vectors on the 
(0, k, 1) reciprocal-space rod, continuum scattering can 
be very well approximated by the "truncated Muller 
ansatz" (TMA) function^ 

aA [1 — cos(ttI)} 



u 2 -A 2 c± -^J 2 sm 2 (Trl) ) (25) 



The TMA is plotted in a thin dotted lined in Figs. 0a- 
c for a direct comparison to our chain-MF/RPA result. 
Conveniently, given the experimental resolution width, 
the two functional forms are almost indistinguishable. 
This fact is illustrated in Fig. || that shows the trans- 
verse continuum obtained in the actual chain-MF/RPA 
calculation for q = (0,0,1) (solid line), along with the 
form |25| (dotted line), both profiles being numerically 
convoluted with the resolution function of the instru- 
ment. Resolution effects taken into account, an almost 
perfect match between the chain-MF/RPA calculation 



for BaCu 2 Si207 and Eqs. (|23| ) and fl25[ ) can be obtained 
for the entire range of energy and momentum transfers 
covered in our experiments by choosing A Ci j_ = 5.0 meV 
and a — 0.17 meV -1 . 
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Just like the transverse-polarized part, the fit function 
for longitudinal scattering is composed of a single-mode 
and a continuum components. The dispersion relation 
and dynamic structure factor for the single-mode contri- 
bution are written as: 



s| M (q,a;) 



7 [I-cos(ttZ)] 
2 



T/tt 



2 wn(q) [uj 

-J 2 sin 2 (nl) + A 2 



"II fa)] + r 2 
A 2 J'(q) 
\J'\ 



(26) 
(27) 



These equations are a generalization of Eqs. (12) and 
(13) in Ref. ^, that allow for a damping of the longitudi- 
nal mode. The adjustable parameter A|| is the the energy 
of the longitudinal mode at the RPA "magic" point. The 
coefficient 7 is an adjustable parameter that determines 
the intensity ratio of longitudinal and transverse excita- 
tions, while A is an overall intensity prefactor used for 
both polarizations (see Eqs. (4), (5), and (10) in Ref. ||). 
In Eq. (|2^) the (^-function is replaced (for positive energy 
transfers) by a Lorentzian profile with a half-width at 
half height of T. 

The longitudinal-polarized excitation continuum was 
modelled using the same truncated Muller-ansatz cross 
section function as previously done for the transverse 
case: 



sj(q,«0 



f3A [1 - cos(ttZ)] 



CJ 2_ 2 H J2sin 2 ((?||) 



-J 2 sin 2 (7r0 (28) 



Note that, unlike in Ref. ||, we use separate relative in- 
tensity prefactors and (pseudo)gap energies for the trans- 
verse and longitudinal continua. By choosing A C) || = 
A Cj j_ and (3 — a one can accurately reproduce the chain- 
MF/RPA result for BaCu2Si2C>7 to within resolution ef- 
fects in the energy and momentum transfer range covered 
in the present study. 



B. Transverse polarization 

As a first step in the quantitative data analysis, the 
two-component model cross section for transverse polar- 
ization was numerically convoluted with the calculated 
spectrometer resolution function and fit to the trans- 
verse components of all scans measured in this work 
(429 total scan points). The relevant parameters of the 
model, including the mass gap A = 2.51(2) meV, the 
inter-chain exchange constants J x = —0.460(7) meV, 
Jy = 0.200(6) meV, 2J 3 = 0.152(7) meV, the in-chain 
exchange parameter J — 24.1 meV, the continuum gap 
A c ,_l = 4.8(2) meV, and the ratio a = 0.20(3) meV -1 of 
single-mode and continuum intensities- were determined 
previously with very good accuracy.El£2l In analyzing the 
present data, only the overall scaling factor was treated as 



an adjustable parameter. A good (x 2 — 2.7) 1-parameter 
global fit to all the measured scans was obtained (heavy 
solid lines in Figs. ||a-^a and Figs. ||a-g. The hatched and 
greyed areas represent the continuum and single-mode 
components, respectively. 

As mentioned in the previous section, our chain- 
MF/RPA theoretical result for BaCu2Si2C>7 corresponds 
to fitting function parameters A Ci j_ = 2A = 5.0 meV and 
a = 0.17 meV -1 , which is in a remarkably good agree- 
ment with previous and current experiments. We con- 
clude that for transverse polarization the chain-MF/RPA 
not only predicts the correct spin wave dispersion rela- 
tion and continuum gap energy, but provides an accurate 
estimate for the intensity of the lower-energy part of the 
continuum. 



C. Longitudinal polarization 

The agreement with theory is not nearly as good in the 
longitudinal polarization channel. In the chain-MF/RPA 
the LM is infinitely sharp and corresponds to T — > in 
Eq. (p6|). The LM's energy and intensity are given by 
A|| = V3A, and 7 « 0.49. Our chain-MF/RPA calcu- 
lation for the longitudinal continuum in BaCu2Si207 is 
very well approximated by Eq. (28) with A c .^ = A c ii 
and f3 = a. Using these values in the model cross section 
convoluted with the resolution function of the spectrom- 
eter, we can simulate the measured scans as expected in 
the chain-MF/RPA model. These simulations are shown 
in solid lines in Figs. |2f3-^b. The dark greyed area repre- 
sents the longitudinal mode, and the hatched area is the 
continuum contribution. It is clear that at all values of 
q_i_ the model fails to reproduce the observed longitudinal 
spectrum. The discrepancy is greatest at energy trans- 
fers below 2 A, where the chain-MF/RPA model predicts 
no scattering except that by the LM. At higher energy 
transfers the agreement between theory and experiment 
becomes progressively better. 

Of course, much better fits to the experimental data 
can be obtained if the central energy Ay, intensity pref- 
actor 7 and intrinsic energy width T of the longitudi- 
nal mode are allowed to vary. The result of fitting this 
"damped LM" model globally to the entire data set for 
longitudinal polarization (358 data) is shown in Figs. ^]b- 
||b in a dotted line, and corresponds to \ 2 — 1.5. The fit 
yields Ay = 2.1(1) meV, 7 = 1.2(2) and T = 1.5(2) meV. 
This analysis confirms the main conclusion of the prelim- 
inary study of Ref. ||: to adequately describe the longi- 
tudinal scattering in BaCu 2 Si2 0y in terms of a "longi- 
tudinal mode" one has to assume a substantial intrinsic 
width, comparable to the mode's central energy and to its 
separation from the continuum threshold. The "longitu- 
dinal mode" can therefore be no longer considered a sepa- 
rate feature, as it is merged with the strong continuum at 
higher energy transfers. The energy separation of single- 
mode and continuum excitations previously observed for 
transverse polarization is absent in the longitudinal chan- 
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nel. It is important to emphasize that the mismatch be- 
tween theory and experiment involves more than simply a 
broadening of the LM. Experimentally one observed con- 
siderably more scattering below 2A energy transfer than 
the LM could provide in the chain-MF/RPA model. As 
a result, the refined value of 7 is almost 4 times larger 
than expected, and the "LM" is almost equal in intensity 
to a transverse spin wave. 

The measured data can, in fact, be reproduced without 
including a single-mode longitudinal component in the 
cross section. This "continuum-only" model corresponds 
to 7 = 0, while A c || and a\\ are the adjustable parame- 
ters. Rather good global fits to 249 data points at k = 
are obtained with A c ,|| = 2.0(1) mcV, a\\ = 0.22(0.01) 
and x 2 = 1-16. Scan simulations based on these param- 
eter values are plotted in heavy solid lines in Fig. [|b,d,f, 
and h, and in a dash-dot line in Fig. ||b. The param- 
eter A c || was fit separately for the constant-Q scans at 
k = —0.5 and k = —1, yielding A c 11 = 1.8(1) meV and 
A c || = 1.5(1) meV, respectively. The results are shown 
in dash-dot lines in Figs. ||b and |||d. The variation of 
A c 11 as a function of qj^ corresponds to the dispersion 
of the longitudinal mode built into the "damped LM" 
model. 



VI. CONCLUDING REMARKS 

Based on the neutron scattering results we can now 
give phenomenological description of the longitudinal 
excitations in weakly interacting quantum spin chains. 
There is no sharp longitudinal mode, but a broad asym- 
metric peak that is inseparable from the continuum at 
higher frequencies. This feature is practically indepen- 
dent of q^, but has a steep dispersion along the chain 
axis. The scattering starts at energies well below 2A, 
and its intensity at low energies is considerably greater 
than predicted by the chain-MF/RPA. 

It appears that the established chain-MF/RPA model 
is at the same time remarkably good in predicting the 
transverse correlations of weakly-coupled chains, and 
sourly inadequate as far as longitudinal fluctuations are 
concerned. Admittedly, one can never entirely dismiss 
the possibility that the disagreement between theory and 
experiment in the latter case may, in fact, be due to some 
intrinsic flaw in the unconventional technique that we 
used for polarization analysis. However, having repeat- 
edly scrutinized the measurement procedure, we were un- 
able to identify any potential sources of systematic er- 
ror that could account for the observed discrepancies 



with theoretical calculations. We thus conclude that 
the discrepancies stem from limitations of the theoretical 
method itself. Among the assumptions and approxima- 
tions associated with the chain-MF/RPA approach, the 
most likely source of errors is the uncontrolled discarding 
of the self-energies in the RPA. The RPA, by definition, 
acts on bare (purely ID) dynamic susceptibilities at par- 
ticular wave vectors. It excludes interactions between 
particles, such as processes that involve a decay of a par- 
ticle with momentum q into a pair of particles with mo- 
menta qi + q2 = q- The contributions of such processes 
to the susceptibility involve ID correlation functions of 
three or more spin operators. For spin chains that are 
intrinsically gapped such processes are expected to be 
suppressed, in which case the RPA will be fully justi- 
fied. We can expect the RPA to be an almost perfect 
description of weakly coupled ladders or Haldane spin 
chains. For weakly coupled 5 = 1/2 chains however the 
mean- field gap A is itself determined by J'. As a result, 
the transverse spectrum in RPA is gapless, regardless of 
J' I J. Hence a longitudinal excitation can always de- 
cay into a pair of transverse-polarized spin waves. The 
RPA fails by excluding this effect. Comparing the re- 
sults of the present study to the ones reported in Ref.|l| 
for KCUF3, the question arises why there is a longitudi- 
nal mode, albeit damped, in the latter material but not 
in BaCu2Si207. The main difference between the two 
materials is the strength of the interchain coupling: in 
BaCu2Si2 07 the ratio of the bandwidths perpendicular 
to the chains and along the chains is 2A/ttJ ~ 0.066, 
whereas it is approximately 0.2 for KCUF3. This may 
suggest that a sufficiently strong dispersion perpendic- 
ular to the chains is necessary in order to stabilize a 
damped longitudinal mode. It would be interesting to 
investigate this issue by determining the damping of the 
longitudinal mode in MF/RPA. 
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